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Abstract 

We consider Mellin-Barnes transform of triangle ladder-like scalar diagram in d — 4 dimen- 
sions. It is shown how the multi-fold MB transform of the momentum integral corresponding 
to an arbitrary number of rungs is reduced to the two-fold MB transform. For this purpose we 
use Belokurov-Usyukina reduction method for four-dimensional scalar integrals in the position 
space. The result is represented in terms of Euler ^-function and its derivatives. We derive 
new formulas for the MB two-fold integration in complex planes of two complex variables. 
We demonstrate that these formulas solve Bethe-Salpeter equation. We comment on further 
applications of the solution to the Bethe-Salpeter equation for the vertices in M — 4 supersym- 
metric Yang-Mills theory. We show that the recursive property of the MB transforms observed 
in the present work for that kind of diagrams has nothing to do with quantum field theory, 
theory of integral transforms, or with theory of polylogarithms in general, but has an origin in 
a simple recursive property for smooth functions which can be shown by using basic methods 
of mathematical analysis. 
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1 Introduction 



In this paper we study Mellin-Barnes (MB) transform of three-point scalar ladder- like integrals in 
d = 4 space-time dimensions. These integrals contribute to Green's functions of four-dimensional 
massless scalar theories and can be represented in terms of UD functions [TJ [5] . As has been proven 
in Refs. [U |31 [S] at the level of graphs, the UD functions transform to themselves under Fourier 
transform. Also, this property can be proven by making a use of the MB transform [BJ. 

As has been shown in Ref. [5], any other scalar three-point Green's function even in the non- 
integer number d = 4 — 2e of dimensions possesses such a property of invariance with respect 
to Fourier transform too. This happens due to the fact that any three-point function can dia- 
grammatically be considered as a net of consequent three-point or four-point integrations in the 
position space since any of those integrations in the position space can be represented in terms of 
MB integral transforms with powers of spacetime intervals in denominators depending on the MB 
transform parameters. Doing consequently momentum integration via the MB transform, we come 
to an expression in which there is no momentum integral but there are integrations in complex 
planes of the MB transform parameters. We can apply arguments of Ref. [5J |BJ to the expressions 
of such a type in order to demonstrate their invariance with respect to the Fourier transformation. 

In general, technique of the MB transform is a powerful method to perform complicate multiloop 
calculations in quantum field theory [7J EJ [9j HO] . For example, the three-point one- loop function 
in the position space in the massless scalar theory is a combination of Appell's hypergeometric 
functions which have appeared as a result of residue calculation via the MB transform. In Ref. 
[6j it has been shown that for the three-point integrations in the position space with the powers of 
space-time intervals in the denominators equal to the integer numbers shifted by some multipliers 
of e, where e is a parameter of dimensional regularization, d = 4 — 2e, these Appell's hypergeometric 
functions can be reduced to a combination of the UD functions. On the other side, the four-point 
function in the position space in the scalar theory in d = 4 dimensions can be represented in terms 
of the UD functions too |12j . The invariance of those four-point ladder integrals under the Fourier 
transform can be proved by means of the MB transform too [12]. All these properties make UD 
functions very attractive for further investigation by making a use of the MB technique. The main 
advantage of this technique that we use in the present paper is that the three-point d-dimensional 
integration is transformed to a power-like form, which can be integrated in further loops as a three- 
point integration with the powers of space-time intervals in denominators which depend on the 
complex variables of the MB transform. 

The first aim of this paper is to investigate further the properties of the UD functions via the MB 
transform. As a by-product, we derive new formulas for the MB two- fold integration in the complex 
planes of two variables of integration. The third aim is to demonstrate compatibility between these 
formulas for the MB transform and the Bethe-Salpeter equation for the infinite sum of triangle-like 
and box-like scalar ladders. 

The fourth motivation for the investigation developed in this paper is to apply the experience 
obtained with studying the Bethe-Salpeter equation for the scalar ladders to calculation of an 
auxiliary vertex or Green function, which in Af = 4 supersymmetric Yang-Mills theory does not 
depend on any scale, ultraviolet or infrared, in the limit of removing dimensional regularization, 
e = 0, in d = 4. That vertex should satisfy the Bethe-Salpeter equation, particular for this Green's 
function. This is Lcc vertex in which an auxiliary field L couples to two (self-adjoint) Faddeev- 
Popov ghost fields c. To investigate that vertex by the Bethe-Salpeter equation has a sense since 
the ghost fields are the only fields of this vertex which stand in the measure of the path integral, 
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and therefore this equation restricts the vertex strongly. 

As has been shown in Refs. |T3l HH HI], in Af = 4 supersymmetric Yang-Mills theory starting 
with this vertex all the other correlators can be found in the four-dimensional limit (when d = 4) 
by using Slavnov-Taylor identity [HI H3 [HI [13 [20l EU [22l [23]. The explicit two-loop result for 
that vertex has been found in Refs. [24l [25] [26] . All the poles in e disappear in all number of loops 
for this vertex, however it cannot be analysed by the methods of conformal field theory since the 
auxiliary field L does not propagate. Nevertheless, the three-point function of dressed mean gluons 
of Refs. [131 d3J [TS] can be fixed by conformal symmetry in analogy to [571 HH 12S1 1301 131] > however, 
we cannot find higher-point correlators of the dressed mean fields from this Green's function by 
using ST identity. 

As it can be seen from Refs. [T31 dH [13], the explicit structure of Lcc correlator includes 
logarithms and the UD functions of ratios of spacetime distances between arguments of the Green's 
functions. As it follows from Ref. [4], the ladder-like diagrams with horizontal gluon lines, that are 
important subclass of all the diagrams contributing to that vertex can be represented in terms of 
the UD functions too for an arbitrary number of the horizontal gluon lines. In the next paper we 
will demonstrate that any of the contributions to the correlator can be represented in terms of the 
UD functions [32]. 

In this paper we consider a particular solution to the Bethe-Salpeter equation, that represents 
an infinite sum of the triangle scalar ladder diagrams in four space-time dimensions in the MB 
representation. This particular solution has already been studied in the case of the box-ladder 
diagrams by Broadhurst and Davydychev in Ref. [31] ■ The authors of [31] calculated a sum of 
those box-ladder diagrams by making a use of the integral representation for the UD functions of 
Ref. [2]. As has been shown in Refs.[IJ[2], the triangle- like and the box- like ladders in four space- 
time dimensions are related by conformal transformation in dual space and their MB transforms 
coincide. The difference of the present paper with the approach of Ref. [34] is that we do the 
calculation via the MB transformation. We derive recursive relations for the MB transforms of 
the momentum integrals corresponding to the UD triangle diagrams by using Belokurov-Usyukina 
trick [33] for the ladder diagrams in the position space, reducing the multi-fold MB integrals to the 
two-fold MB integrals. 

2 Loop reduction in d = 4 dimensions 

In this section we describe a trick with a help of which the triangle ladder diagram of n loops is 
reduced to the triangle ladder diagram of n — 1 loops. This trick has been published for the first 
time in paper [33 a . The effect of the loop reduction has earlier been discovered in Ref. [35] for 
a propagator-type diagram of a particular topology in a special limit for the indices of diagram. 
Note, that in Ref. [33] the derivation has been performed without that special limit. Later, the key 
intermediate point of the trick has been published in Ref. [35] . Since all of that together has never 
been published in one article in detail, we collect all the results and intermediate steps we need for 
deriving our formulas for the MB transform of the triangle-ladder diagram. 

The basic formulas we used to derive these diagrammatic relations are the chain integration for 
an any <x\ and «2, except for d/2 + n, n G N, 
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Figure 1: Basic relations 
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Figure 3: Derivation of the basic relation. Part II. Lines without any index are the lines with index 1 
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Figure 4: Transformation of the first diagram in Fig. (1). Lines without any index are the lines with index 1 



A new d-dimensional measure, introduced in ref. [24], 

Dx = Ti^id d x 

is used to avoid powers of it on the r.h.s. of the diagrams. 

The bold dots correspond to the internal vertices of the diagrams, and this measure corresponds 
to the bold dots. We make an integration over coordinates of the internal vertices. In the first line 
of fig. © we present the formula which we want to derive. This formula is present in paper [36 , 
however with a bit different indices of the diagrams. Later we will see how the r.h.s. of the first 
line of fig. © can be put in equivalence to fig. (4) of [35]. The second and the third lines of fig. © 
are consequences of the first line. How exactly the second and the third lines can be derived from 
the first line we explain further in the present paper. Usually, it can be done by inserting points of 
integration into the lines of the graphs. The third line is published in Ref. [33]. The second line 
has never been published. The condition for the e-terms in the indices is 

£i+ e 2 + £3 = 0. 

In fig.© and in fig.©, the latter is a continuation of fig.©, we have shown how the first line 
of fig.© can be derived. Actually, the only mathematics that has been used here is the uniqueness 



6 




r(i - £l )r(2 + £ i) r(i - £3 )r(i - £2 )r(- £l ) 



r(i + ei)r(- £l ) r(i + £3 )r(i + £2 )r(2 + £l ) 



2 + El 




1-ei J 



r(i- £l )r(2 + £l ) 
r(i + £l )r(- £l ) 



r(l-ei)r(2 + ei) r(-ei)r(l+ei) 



r(i + ei)r(-ei) r(2 + £l )r(i- £l ) 

l + ei 



1 + e 3 J 





r(2- £2 )r(i + £2 ) r(i + £l )r(2 + £2 )r(i + £3 ) 
r( £2 )r(i - £2 ) r(i - £l )r(- £2 )r(i - £3 ) 



+ £2 



1 + £3 J 



r(2 - £2 )r(i + £2 ) r(i + £l )r(2 + £2 )r(i + £3 ) r( £2 )r(i - £2 ) r(- £2 )r(i + £2 ) 



r( £2 )r(i - £2 ) r(i - £2 )r(- £2 )r(i - £3 ) r(2 - £2 )r(i + £2 ) r(2 + £2 )r(i - £2 ) 



i + £2 



l + £l 




r(i + £3 )r(2- £3 ) 
r(i - £3 )r( £3 ) 



i + £2 



i + £3 



1 + 1 



i + £2 



i + £2 



l + £ 3 

1 +El 



l+e 3 



r(i + £3 )r(2- £3 ) r( £3 )r(i + £l )r(i + £2 ) _ \ + £ ; 



l+£l 



l + £2 



1 +El 



r(i - £3 )r( £3 ) r(2 - £3 )r(i - £l )r(i - £2 ) 



i + £2 



i + £2 



l+e 3 

1 + El 



1+E3 J~ l 
1 + £2 



Figure 5: Transformation of the second diagram in Fig. (1) 




Fi gure 6; Fourier transform of the diagrams in Fig. (1). The r.h.s. corresponds to the momentum integrals with 
indicated indices. 



star-triangle relation and the integration by parts, both the relations are well-known in the scalar 
field theories [37l |38l [39] (for a short review, Ref. [40]). In all the figures from ([I} till (|9]) we use 
for factor J the definition 

r(i- £l )r(i-g 2 )r(i-£ 3 ) 
' r(i + Ei)r(i + e 2 )r(i + e 3 )' 

taken from Ref. [3"3"] . 

In fig. (j4]) we show how the first diagram on the r.h.s. of the first line of fig. ([1]) transforms 
to the first diagram on the r.h.s. of the second line of figure ([I}. Again, we converted non-unique 
vertices to unique vertices by inserting the point of integration into one of the propagators attached 
to the vertex. Then, the star-triangle relation has been applied in a direct way or in a triangle-star 
way, until the diagram is converted to a diagram with desirable values of indices. The line without 
any index means that the index of the line is 1. We put the exact value of this index on the 
corresponding propagators of fig. ([!} only. In the rest of figures we omit the index 1 just not to 
overload the diagram with indices. 

In fig. ([5]) the second diagram of the first line of fig. ([TJ is transformed to the second diagram of 
the second line of fig. ([IJ . Again, nothing else but the creating of unique vertices and triangles by 
putting new points of integration into the propagators has been used to obtain the desired indices. 
Thus, we have proven the first line and the second line of fig. (Q}. 

The third line of fig. ([1]) can be obtained from the second line of fig. ([1]) by simple integration. 
Indeed, we eliminate the leftmost propagator on the l.h.s. of the second line by convoluting its 
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Fi gure 7' Fourier transform of the diagrams in Fig. (1) 
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Figure 8: Formula (25) of Ref. Q] 
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To derive formula © we need to replace in Eq. (JXJ) each factor in the integrand with its integral 
Fourier transform. The formulas for the Fourier transform of the factor like those can be found in 
Ref. [40]. Indeed, since 

J dpe { P 2 r~ r(a) UV 



1 _--d/2,-q r(d/2-a) 
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As a consequence, we obtain in case of ot\ + a 2 = d 



2.1 Description of Fourier Transform for the diagrams 

In fig. (j6j) the Fourier transform is made for all the diagrams that appear in the first line of fig. ([T]). 
The procedure of Fourier transformation is done as follows. We replace each factor in the integrand 
of the position space representation with the integral Fourier transform of the corresponding factor 
in the momentum space representation. Space-time coordinates appear in the exponentials of the 
Fourier transforms. Integrating over the coordinates of internal vertices we create Dirac (^-functions, 
corresponding to the momentum conservation in each vertex of integration in the position space 
(internal vertex). The momentum integrals over loop momenta will be the Fourier transforms of 
the integrals in the position space. For example, the described procedure in a particular case of the 
one-rung integral in the position space would result in 

[31] I ^ dAz [2y][ly][lz][yz][2z] = (2^ / ******* 5 ^ +») X 

x e <P2* 2 gipia.1 j P3 x 3 C (2) (p2 ) p 2 p 2 y 

We assume the notation of Ref. [24 , where [Ny] = (xjy — y) 2 and analogously for [Nz] and 
[yz], that is, N — 1,2,3 stands for xjq — x\,x 2 ,x^, respectively, which are the external points of 
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the triangle-like ladder diagram. In Refs. [TJ [2] C^(pl,pl,pl) is the definition for the result of 
momentum integrals for n-rung triangle ladder diagram in the momentum space representation. 
Just to make clear the definition for the operation of Fourier transformation used in fig. (|6|), we 
provide an example of the relation 
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The transformation from the position space to the momentum space is necessary to normalize on the 
results for the MB transform of the UD functions done in [I] [2J for the momentum space integrals. 

Fig. ([7]) is fig. ([6]) with signs changed for all e-terms. After this change of signs, the Fourier 
transform of the first line of fig. ([IJ takes the form depicted in fig. ©. In that diagrammatic 
relation we recognize Eq. (25) of Ref. [T]. 

3 Recursive formula for MB transform of the UD functions 

We start this section with writing the definition of the momentum integral for the diagram in the 
r.h.s. of Fig. (|5J). To define the integral measure in the momentum space we use the notation 



Dk = n-*d a k. 



(3) 



This is done to avoid powers of ir in formulas for the momentum integrals. For example, fig. (|8j) 
takes the form shown in fig. (|9j). 
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Figure 9: Fig. (O after measure redefinition 10 

We take a formula from paper [I] for the representation of the three-point momentum integral 

1 



J(ui,V2,u 3 ) = J Dk 
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in terms of the MB transform 
r (-z 2 -V2-V3 + d/2) r (-z 3 - i/! - 1/3 + d/2) r (z 2 + z 3 + r (Ei/< - d/2 + z 3 + z 2 )} = (4) 

1 



(pl)^- d/2 j c 

where we have introduced the definition of function D, 

D ( Z2 ,z 3 )< v v 1 _ r (-z 2 ) r (-£3) r (-z 2 - ^ - ^ + d/2) r (-z 3 - ^ - ^ + d/2) 
[1 ' 2 ' 3J n,r(y.) 

r (z 2 + z 3 + z/ 3 ) r (Si/j - d/2 + z 3 + z 2 ) 

r(d-s^) ■ (5j 

Here and in the rest of the paper we use the notation 

2 2 

= pi = n L 
x — 2 ' y — 2 ' 

P3 Pi 

where the d-dimensional momenta pi,P2,P3 satisfy the conservation law p\ + p 2 + p 3 = and are 
related to d-dimensional momenta q±, q 2l q 3 by a parametrization 

Pi = 93 - 92, 

P2=Ql- 93; 

Pa =92-91. 

As a consequence of this definition, pi appears to be a momentum that enters the one-loop triangle 
diagram in the vertex of the triangle which is opposite to the line with index v\. 



dz 2 dz 3 x»y*>D('* z *)[v 1 ,V2,va] 



3.1 Description of notation in the MB representation 

We absorb into the definition of the MB transform Z)( Z2 > Z3 ) [ui, v 2 , v 3 ] of this three-point integral 
all the factors except for a power of the square of the external momentum p 3 . For d = 4 in the 
denominator we have a sum of the indices minus two for the power of p 3 . 

We do not write the powers of i, supposing that we work in Euclidean space and the corre- 
sponding power of i can be recovered back after Wick rotation. 

The contour of integration C passes a bit on the left of the imaginary axis, separates left and 
right poles and should be closed to the left infinity or to the right infinity. We choose to close 
the contour of integration in the complex plane to the right infinity. It could be closed to the left 
infinity too but it makes more complicate to take the residues into account since the residues in 
variables z 2 and z 3 are mixed in that case. Whether we have to close the contour to the right 
infinity or to the left infinity, the result should be the same function. We omit the factor l/2ni 
that accompanies each integration over MB transform parameter. The inverse factor is generated 
in front of the residues. 

It is known that this representation of the three-point integral can be derived by applying the 
two-fold MB transform to the integral over Feynman parameters, with producing the Euler beta 
functions |7j after integrating over these parameters. There is a difference with the representation 
used in Ref. [5] . The form of Ref. [B] can be recovered from (U) by a cyclic redefinition. In the rest 
of the paper, we use representation 
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3.2 Multi-fold MB transforms of the UD functions 



To develop the recursive equations for the MB transforms of the UD functions, we recover the 
definition of the MB transforms for the triangle ladder diagrams in d = 4 spacetime dimensions 

C^(pI P Ip 2 3 ) = T-L;& n) (x,y) = I dz 2 dz 3 x^y^ M^\z 2l z 3 ), 

\P 3 ) \Pz) JC 

where we used the notation Ai^ (z 2 , z 3 ) for the MB transform of the UD function <&M (x,y) of 
two arguments. The definition of incoming momenta for the triangle ladder diagrams is like in 
Ref . [H [21 [5] . We include in the definition of Ai^(z 2 , Z3) all the factors that can appear in front of 
the MB transformation except for the power of the momentum p 2 . 

The two-fold MB transform is known for the first UD function only. For a UD function of a 
higher order, the result is given in terms of the multi-fold MB transform [lj[2]. In the next section 
we reduce the multi-fold MB transform to the two-fold MB transform by making a use of the loop 
reduction trick described in the previous section. 

The iterative integral relation for the triangle ladder diagram given in Refs. [U El H] , 

C<»V P 2 P 2 ) - f d\ C^((Pi+r n f,(P2^r n r, P l) 

results in an integral recursive relation for the MB transforms of the UD functions 

$ („+l) (p\ PT] 2 / d 4 r $ („) ( (Pi+r) 2 fe-r) 2 \ 1 

\pVpIJ V V Pi ' Pi /(P!+r) 2 (R-r)V 



P i 1 *r 6 dz 2 dz 3 M^(z 2l z 3) i^ +r)2 y 2 f^-^y 3 ' 



c 

d A r & dz 2 dz 3 M {n \z 2 ,z 3 ) 



Pi ) V P% ) (Pi + r) 2 (p 2 - r) 2 r 2 
1 1 



((pi +r) 2 y- Z2 ((p 2 -r) 2 ) l -^r 2 (p 2 ) 2 ^- 1 

TT 2 <j) dz 2 dz 3 M [n \z 2 ,Z 3 )J(l ~Z 3 ,l- Z 2 , 1) ^2)^3-1 = 

7T 2 <j) dz 2 dz 3 du dv x u y v M {n \z 2 ,z 3 )D {u > v) [l - z 3 ,l -z 2 ,l], 

where D^ u ' v > [1 — Z3, 1 — z%, 1] is defined above as the MB transform of the corresponding three-point 
integral, according to description done in the previous subsection 0. 

Thus, one can derive the first recursive formula for the MB transforms, which relates the trans- 
forms of two neighbour UD functions [TJ [2] 

M (n+1) (u,v) = ir 2 / dz 2 dz 3 M {n) (z 2 ,z 3 ) D("^[1-z 3 ,1-z 2 ,1]. 



The function ^"'^[l - z 3 ,l - z 2 , 1] is a nontrivial combination of the Euler gamma functions 
in numerator and denominator, 

[i - *, 1 - i] = ri V ( ,i ri Z~t7' Z ^n + X V) r r r d + « + • 
r (1 + z 2 + z 3 ) r (1 - z 2 ) r (i - z 3 ) 

This formula is written from definition ((4]). 

1 Factor it 2 appears due to usage of the original measure d 4 r in the momentum space. 
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4 Reduction of the multi-fold MB transforms 



The formulas of the previous subsection are in some sense consequence of the ladder-like topology 
of the diagram. We did not make any integration in the complex planes of MB parameters. 

Consider the diagram on the l.h.s. of fig. ©• We can repeat the trick of the previous section 
and obtain the following MB representation for it 
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_ ' {Pi+rf \ 22 f ( P 2-r) 2 \ Z3 g(*2,*3)[i + g2) i + ei) i + g3 ] 

Pi J f^ 3 \ pi ) \ pi J [{ p 1+r) 2 ] l + e 2[{p2 _ r) 2 ] l +Sl[r 2 ] l + e 3 

1 I 1 D(^»)[l + E2 ,l + ei,l + E3] 



Pi Jc ^ dZ2dZ3 [P2\ Z2+Z3 [(Pi + r) 2 Y+^-^[{p 2 - r )2]l+e 1 - 23 [ r 2]l+6 3 
4 / dz 2 dz 3 2 1 J(l + £i -23,1 + £2 - Z 2 ,l + £ 3 )£> (Z2 ' Z3) [1 + £2,1 + El, 1 + £ 3 ] = 

Ps Jc \Pi\ Z2+ 3 

-^-r l dz 2 dz 3 du dv x u y v D {u ' v) [1 + £i - z 3 , 1 + £ 2 - z 2 , 1 + £3] x 

m) Jc 

£)U2, Z 3)[ 1+e2)1+£l)1+£3 ] (6) 

This formula is derived in analogy with the MB transforms of the UD functions in the previous 
section. In detail, the procedure looks like follows. First, we calculate the MB transform of the 
leftmost triangle integral, this yields another triangle integral with indices depending on the complex 
variables of the previous MB transform. This procedure will be used in all the constructions below. 

The general strategy of the present investigation consists in the expanding of the r.h.s. in terms 
of Si. The coefficients of this expansion will contain ip functions and its derivatives. In the r.h.s. on 
fig. © we cannot put all the values of e, immediately equal to zero without expanding these r.h.s. 
in terms of e% and observing that poles in e% disappear. Thus, instead of a Laurent series we obtain 
a Taylor series. Comparing the coefficients in front of the different products of powers of 64 we can 
derive an infinite number of new relations for the two-fold MB integrals over the complex variables 
z 2 and z 3 . 



4.1 Description of the momentum integral 

Notation to be used 

(£1, £2, £3) 

corresponds to the momentum integral with incoming momenta Pi,p 2 ,p 3 as they are depicted on 
the l.h.s. of Fig.®, with the integral measure defined in (j3|), but with n loopfH. The £-terms 
£1, £2, £3 appear in the indices of lines for the first eight propagators on the left side of the diagram, 
namely, in the indices 1 + £1, 1 + £2, 1 + £1 on the upper side of the diagram, 1 + £2, 1 + £1, 1 + £2 
on the lower side of the diagram, and 1 + £3 on the first two rungs. The rest of lines have indices 
equal to 1. These positions and values of the indices are indicated in the third line of Fig. ([J) (in 
the position space). 

2 with amputated external legs 
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4.2 Two-fold MB transform for the two-rung ladder 

As a first step to our formulas, we reproduce Eq.(25) of Ref. 1 . This is already done in Section 
2, however we make it again just to introduce the notation that will be used for the higher rung 
diagrams. 

We consider a diagrammatic relation that can be obtained by integral convolution of the leftmost 
external point of the diagrams in the first line of fig. ([TJ with the line that has index 2 — £3. In such 
a way the Dirac 5-function is produced which eliminates one of the integrations and the leftmost 
point on the l.h.s. is converted to the external vertex, that is, it is not a vertex of integration longer. 
The formula which can be used for this purpose is identity ([2]). 

Then, we take the Fourier transform of each one of the diagrams in the diagrammatic relation 
obtained in that way. By keeping all the factors that appear after making the Fourier transform on 
both the parts of this new diagrammatic equation, we come to 



J 



3 r(i 



■£3) 



1 



J 



J r(l+e 3 ) 



{— £ lj — e 2, —£3) — 
■7(1,1,1 + 61) + 



J r(i + e 3 ) 



1 



■^1 /.^i -^(pj)l+es J ( 1 ' 1 ' 1 £3 ) + 



eie 2 r(l-e 3 ) (plY^iPi 

J r(i + e 3 ) 



1 



J(l,l,l + e 2 ) 



After a simple algebra we derive an equation 

J- 1 



( — £ 1> ~ e 2, —£3) 



£2S 3 {pD 62 



J(l,l,l + £l) 



1 



£\£2 



-J(l,l,l-£ 3 ) 



£i£ 3 {PiY 



-J(l,l,l + £ 2 ) 



from which by changing the signs of all the values £j we obtain 



j(£i,e 2 ,e 3 ) 



J 



1 



1 



£2£ 3 (p|) 



-J(l,l,l-£l) 



£l£2 



-J(l,l,l + £ 3 ) 



r—J(l,l,l-£2) 



This is exactly Eq.(25) of Ref.Q]. This equation has been derived diagrammatically in the previous 
section. In that case the Fourier transform of the first line of fig. ([I} has been done. This approach 
to reproduce Eq.(25) of Ref. 1 can be considered as a cross-check of the procedure used in the 
previous section to produce fig. ©. 

At this moment we start to use the matter developed in the previous section. First of all, we 
rewrite the r.h.s. of Eq.([7]) in the MB transformed representation 



(£1, £2, £3) 



J 



(plY-es 



(pIY 2 



e 2 £ 3 (pi) 1 "^ 



1 



£i£2 (p§) 



du dv x u t^D^pL+ea] + 



(P? 



2\£l 



£i£ 3 (Pi) 1 " 62 



du dv x u y v D^[l-e 1 ]+ 



du dv x u y v D {u > v \l - e 2 ] 



(8) 
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where we have introduced a notation 



D (u,v) \i + v ] =£>(«■*) [1,1,1 + ^]. 



Eq.© can be written as 



L (ei ' £2 ' £3) = w (iY i du dv xU yVDiu,v)[1 ~ £l]+ 

— I dudvx u y v D^[l+e 3 ] + — (^] <f du dv x u y v D {u ^[l - s 2 ] = 

T -l £ 2 Jc £ 1 £ 3 \P3/ JC 

du dv x u y v+e2 D {u ' v) [l-e 1 } + 



J 



£ 2 £ 3 JC 



1 



J 



du dv x u y°D ( - u ' v) [l + £3] + — / du dv x u+£l y v D M [1 - e 2 ] 
£ 1 £ 2 Jc £ l £ 3 Jc 

£)(«.«-ei)[i_ ei ] [1 + £3] £)("- e i' u )[l-e 2 ]' 



bi) 2 



du x" y 11 



£263 



£l£2 



£l£3 



(P 2 ) 2 JC 



du dv x u y v M^' v) [ £l ,e 2 ,e 3 }. 



(9) 



Here we shift the variable of integration in the complex plane. That means, the contour of in- 
tegration C still passes between the left and the right poles. Positions of the poles in the plane 
of two complex variables will be changed with that trick but their nature (left or right) remains 
unchanged. 

On the other hand, we have formula ([6]), from which we obtain 



j> dz 2 dz 3 [1 + £1 - z 3 , 1 + £2 - z 2 , 1 + £ 3 ]£> (22 ' 23) [1 + £ 2 ,1+ £1, 1 + £3] 



J 



D (U,V-S 2 ) []_ _ £l j £>(".-») [l + £ 3 ] £)(«-£!,!)) ^ _ £2 ] 



£2£3 



£l£2 



£l£3 



£ i) £2, £3] 



(10) 



This formula is valid for any u, v and presents by itself a nontrivial result which can be used in 
practical applications of the MB integration. This is the two-fold MB transform of two-loop integral, 
and it is valid not only in the limit of all £j — > but for any £j . Expanding this formula in terms of 
£i we obtain infinite number of new relations. The explicit form of M 2 U [£1 , £2, £3] is calculated in 
the next section in the limit of vanishing £j, the result is Eq. (jTHJ) , 



£ dz 2 dz 3 D^[l-z 3 ,l-Z2,l]D^^[l,l,l] 



T 2 (l + u + v )T 2 {~u)T 2 {-v) x 



+ _ (^(-^ _ ^ ( _ M ))2 _ (r(1 _ £)r(1 + £)) 



(2) 

6 = 



4.3 Two-fold MB transform for the three-rung ladder 

We can derive a new formula for the two-fold MB transformation going to higher loops in the ladder 
diagrams. We start to work with integral J 3 (ei, £2, £3) which appears on the l.h.s. in the third line 
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of Fig.([T|) and take the Fourier transform of each one of the diagrams of that line. The Fourier 
transform of the l.h.s. can be written as 

f r{i-e 1 ) \ 3 f r{i-e 2 ) \ 3 f r{i-e 3 ) \ 2 1 if, e , 

and the whole diagrammatic relation becomes as 



J 



3 r(i + £3 ) 1 



r(l-e 8 )p?P§ 
J 3 r(l + e 3 ) 1 

£ie 2 r(i-e 3 ) iplY-^iplY-^pl 



{— £l, — £2, — £3) — J 



2 1 r(i + ea ) 



{si, -£2, -£3) + 



£ 2 £ 3 r(i - £ 3 ) (AY'^piivD 1 -^ 
J 2 r(i + E3 ) 1 



£ie 3 r(i- e 3 )p?(^) 1 - Ea (^) 1 - 



(-£1) + 
(-£2). 



We kept all the factors that appear after making the Fourier transform on both the parts of the 
diagrammatic equation. Here we use a brief notation 



(-£l) = / (~£l,0,£i), / (-£2) = / (0,-£ 2 ,£ 2 ) 



After a little algebra we obtain an expansion of the three-loop integral in terms of the two-loop 
integrals, 



J 



(— £1, —£2, —£3) — J 1 
-£2, -£3) 



1 



1 



(pD-^ipl) 1 -* 2 £2£3 

1 1 



{pl)- ei {pl)- e2 Pl £i£2 k 
By changing the signs of all values we obtain an equation 

1 



{pl)- e2 {plY~ ei £i£3 



{-ex)- 
(-£2) 



J {£\,£2, £3) = J 



J- 



1 



{pi) ei (pl) £2 Pl £i£2 



{p\) ei {pi) 1+£2 £2£3 

1 1_ 

(£l,£2,£ 3 j+ (p 2 )£2b2) i +£l £i£3 



(£l)- 
(£2) 



The previous equation can be re-written for the MB transforms, and we obtain the result 



/ (£i,£2,£3) = J 

J3 

-1 1 



1 



1 



£ 2 £ 3 {pi) e Hpl) 3+£2 Jc 



du dv x u y v M^ u ' v \e 1 )+ 



J- 



1 



£i£2 (Pi) ei (p|) £2 (p§) 3 Jc 
1 1 



£i£ 3 (Pl) e2 {pl) 3+ei Jc 



du dv x u y v M 2 >,v \e 1 ,e 2 ,e 3 )+ 
du dv x u y v M 2 iu ' v) {e 2 ) = 



J 



1 (pI^ £i 



1 



_£ 2 £3 \Pl) 3 ) 3 - £3 Jc 
\ £ i / „2\ £ 2 



du dv x u y v M^ u ' v \e 1 )+ 



J' 1 — (4) 1 (4) 2 / du dv * u y v Mi u ' v \s 1 ,e 2 ,s 3 ) 

£i£ 2 \PtJ \pi) {pi) 3 ea J c 
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1 p 3 



du dv x u y v M^ v) {e 2 ) 



2\3-£ 3 



£i£3 \pIJ (p 2 3 ) 3 - £3 Jc 
— I dudv x u - £l y v M^ v \ex) + J" 1 — / du dv x u ~ £l y v ~ e * m\ u ' v) {e u e 2 , e 3 ) 

1 



£l £ 3 JC 



du dv x u y v - E2 M^ v) (e 2 ) 



+- 



£l£3 



-M. 



e 2 e 3 

(u,V+S2) 



eie 2 



(pD 3 - 



du dv x u y v Mp u '{e 1 ,e 2 ,e 3 ). 



Here we have shifted the variable of integration in the complex plane. That means, the contour of 
integration C still passes between the left and the right poles. Positions of the poles in the planes 
of two complex variables will be changed but their nature (to belong to the set of left poles or to 
the set of the right poles) cannot be changed with such a trick. 
On the other side, in complete analogy with Eq.Q we obtain 



j(ei,£2,e 3 ) 



(pi) 2 



Dr <p dz 2 dz 3 
1 



(pi + r) 2 \ / (p 2 - r) 



Pi 



2 \ z 3 



pi 



M 2 (Z2 ' Z3) ( £1 , £2 , £3 ) 



(pi) 2 



[(p 1+r )2]l+e 1 [( ? , 2 _ r )2]l+e 2r 2 

Af< z - Z3) ( £l , £2 , £3 ) 



Dr i dz 2 dz 3 ^ Z2+Z3 ^ + r yY+ei-z 2 [(p 2 - r )2]l+e 2 - Z3r 2 



dz 2 dz 3 - „ 

C [pf] 2 +*=+*3 



J(l + e 2 - z 3 , 1 + £ i - z 2 , l)M [ 2 Z2 ' Z3) (e u e 2 , e 3 ) 



(Pi)" - JC 
Thus, we derive that 



^ y dz 2 dz 3 ducfe ajV0 (,,>tO [l + £a -«3.1 + ei - 22, l]M^ 2 ' Za) ( £l , e 2 , e 3 ). (11) 



M^Iei.ea.ea] = £ dz 2 dz 3 Z?^[l + £2 - z 3 , 1 + £l - z 2 , l]M 2 (z2 ' Z3) ( £l , £2 , £3 ) 



J 



£2£3 



£l£2 



£l£3 



(12) 



This formula is valid for any u,v and for any £\,e 2 . Expanding this formula in terms of £j we can 
obtain inhnitc number of new integral relations. In the next section we calculate the limit Si — > 
of M^% x ,e 2 ,e 3 \. 



4.4 Two-fold MB transform for the four-rung ladder 

In the next paragraphs of this section we consider the case of the four-loop momentum triangle 
ladder diagram. It looks like the diagram on the l.h.s. of the third line of Fig.(fT]) but with one more 
rung. To get that graphical representation, the diagrams in the third line of fig. ((T]) are integrated 
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with three more propagators, the index of each propagator is equal to 1. The Fourier transform of 
the l.h.s. of the diagrammatic relation obtained in such a way contains the integral / 4 (ei, £2 ; £3)) 



r(l- £l )\ 3 fT(l-e 2 )\ 3 fT(l-e 3 )\ 2 1 1 



T(l+e 1 )J \T(l + s 2 )J \T(l + e 3 )J vlvlh 
and we derive the identity 



. a r(l + e 3 ) 1 



J 3 r(i + e 3 ) 1 
ei£ 2 r(i - £ 3 ) pIpIpI 



{— El) ~£2: — £3) — >/ 2 



y"(-£i,-£2,-£3) + 



£2£3r(l-£ 3 )p?p|(^) 1 -^ 

J 2 r(l+e 3 ) 1 



£ 1 £ 3 r(l-£ 3 )p 2 ^(p2M- £ 



i-ei) + 



After a little algebra we obtain an expansion of the four-loop integral in terms of the three-loop 
integrals, 

1 1 



y"(-£i,-£2,-£3) = J' 



,2U-e 2 



£2£3 73 



(-ei) + 3 ' 



P 3 £ 1 £ 3 73 
1 



(— £i ; ~£2, —£3) 
1 



(Pi) 1 -" 1 £i£3 J 3 

from which by changing the signs of all the values £j we obtain an equation 



J (£i,£2,£3) = J 



(pD i+s2 £ 2 £ 3 y 3 



(£1) 



J" 



1 



P3 £ 1 £ 3 J 3 
1 



(£i,£2,£3)+ 
1 



£l£3 J 3 



(S 2 ) 



This equation can be re-written for the MB transforms 

1 1 



y"(£i,£2,£3) = J 



£ 2 £ 3 (Pi) 1+£2 (PI) 3+£1 



dv x u y v M^ v \e 1 )+ 



J 



— / dM dw x " y v M 3 iu,v) (e 1 ,e 2 ,e 3 )+ 

du dv x u y v M^ v \e 2 ) = 



£i£ 2 (p§) 
1 ] 



J 



(p§) 4 " £3 



dit cfo a;" y l 



eie 3 {pl)^{plY+^-J c 
J 



£2£3 



M^fa) + ^—A4 u ' v \ £l ,s 2: e 3 ) + ■j—M 3 u,v \e 2 ) 



£l£2 



l£3 



(pi) 



y dw dv x" 2/"Mi tl ' 1 ' ) [£i,£2,£3] 



On the other hand, in analogy to Eq.©, 



(p§)3-e. 



Dr <p dz 2 dz 3 



(jh±_r)'' 
Pi 



2\ 2 2 



(P2 - r ) ; 

pi 



2\^3 ^(z 2 ,Z 3 ) 



(£l, £2, £3) 
^(£l, £ 2 , £3) 



(pi + r) 2 (p 2 - r) 2 r 2 
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1 



n I d d - : M^ Z3) (e u e 2 ,e 3 ) 

Jc IpI\ Z2+Z3 [(Pi + rfY-H(P2 - r)2]i-a r 2 

± T -J(l- Z3 ,l- Z2 ,l)Mt' Z3 \e 1 ,e 2 ,e 3 )^ 



(R 



' ' dz 2 dz 3 dudv x u y v D^[l-z 3 ,l-z 2 ,l}M^ 2 ' Z3 \e u e 2 ,e 3 ). 



We obtain 



J 



£ dz 2 dz 3 D^[l - z 3 ,l - z 2 ,l]M^ Z3 \s u e 2 ,s 3 ) = 

I,., „,1 . 1 {,. .,1 . C, .,1 . 



-2^3 £1^2 £l£3 



(13) 



, .... . . ., . ._ ... . . ... =M { r\e^e 2 M- (14) 

£ 2 £ 3 £l£2 £l£3 

New integral relations can be derived by expanding this formula in terms of £j, and this formula is 
valid for any u, v. 

5 One-fold MB transforms 

In this section we collect the useful formulas for the one-fold MB transforms, 



/ dz x z r 2 (-z)T 2 (l + z) = —^, 
Jc l ~ x 

j> dz x z {T 2 (-z)T{z + l)T* (z)} = - Q In 2 x + 2((2)^ + ln(l - x) \ax + Li 2 (x), 

<j> dzx z T 2 {-z)T 2 {l + z)i}j(l + z) = 



(1-x) 



[-ip(l) lnx + ln(l - x) lax- C(2) + Li 2 (a;)] . 



It is possible to work with the MB transforms by making a use of such tricks like derivation with 
respect to parameter and integration by parts in the complex plane. Here we demonstrate several 
examples. First of all, after the derivation with respect to x of 



we obtain 



T ^- = £ dzx z r(-z)T(l + z) 



= <j> dz x z T{\- z)T{\ + z). 



(l + x) 



The result can be checked directly by the counting of the residues. Now, let us check that the 
integral of the total derivative is equal to zero. Indeed, we have 

0=/ dz -^-{x z T(-z)T(l + z)} =lnx I dz x z T(-z)T(l + z) + 
Jc dz J c 
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dz x z T(-z)T{l + z)(-ip(-z) + tp(l + z)) = 

Jc 

in.,- / dz ^ r(-z)r(i + z) + r dz x z T (_ z j T ( 1 + z ^ 1 + z }_^ 1 _ z }y ( 15 ) 



1 + X 



c 



c 



We calculate all the integrals on the r.h.s. explicitly and show that their sum is zero. The contri- 
bution of residues leads to 

dz a; 2 r(-z)r(l + z)tp(l + z) = - m/>(2) + x 2 i/»(3) + . . . , 

dz x z T(-z)T(l + z)i){l -z) = ip(l) - xip(l) + x 2 tp(2) - x 3 ip{3) + . . . - In at + Y~> 

I dzx^- Z ^ 1 + z) 
Jc z 

Substituting these results in (|15|) we reproduce zero in the l.h.s. of it. 



= In a; — ln(l + x). 



6 Explicit results of the two-fold MB transform 

In order to use formulas obtained for the two-fold integration explicitly, we need to calculate the 
r.h.s. of the equations in the limit of vanishing £j. First, we consider the result for J 2 (ei, £ 2 , £3). 
There are two ways to derive the result in this limit for the r.h.s. 



6.1 Belokurov-Usyukina expansion for the higher rung ladder 

According to Eqs. © and ([3]), we have to find the limit 

£2— >0,£l— >0 

r D (u,v- E2 ) m 1 flMfi + fjl D(*-ei>«)[i -£ 2 ] 
hm J H H - — 

£2^0, £1^0 £ 2 £ 3 £l&2 £1^3 

1 r(-u)r(-v)r(ei + e 2 - -»)r( £l + £2 - i;)r 2 (1 - gl - s 2 + u + g) 
ei£ 2 r(i-e 3 )r(i + £ 3 ) 

1 r(-u)r(e 2 - «)r(ei - u)r(ei + e 2 - ^)r 2 (i - gl - g 2 + u + v) 
~£2(ei + £2) r(i- £l )r(i + £l ) 

1 r(ei - u)r(-«)r(ei + £ 2 - u)T(e 2 - v)T 2 (l ~e 1 -e 2 + u + v) 



lim 

£2->0,£l^0 



£i(£i + £2) 

T 2 {l + u + v) 



lim 

£2— 5-0, £1^0 



r(i - e 2 )r(l + e 2 ) 

1 r(-it)r(-?;)r(ei + £ 2 - it)r(ei + £ 2 - v) 

£i£ 2 r(i-e 3 )r(i + £ 3 ) 

1 r(-u)r(e2-^)r(ei-M)r(£i + £ 2 -^) 
^£ 2 (£i + £ 2 ) r(i -£i)r(i + £i) 

1 r(ei - «)r(-«)r(ei + e 2 - u)r(s 2 - «) 



£l(£l + £ 2 ) 



r(i-£ 2 )r(i + £ 2 ) 



21 



After a simple calculus, the limit is 

lim M^ u ' v} (e u e 2 , e 3 ) =T 2 {1+ u + v)T 2 {~u)T 2 (-v) x 

£1^0, £2^0 

i W(-v) + ^(- U )) - i - V(-u)) 2 - 2 (T(i - e)r(i + £ ))fj 



(16) 



The last term is a value of Riemann zeta function multiplied by the first UD function. This result 
can be represented in other forms. For example, by integrating by parts in the complex u and v 
planes we reproduce result (|18l) of the next subsection. However, this form is important since it 
stands in the right hand side of Eq. (fT2"j) in the limit of vanishing 

6.2 Integration by parts in the complex plane 

Also, it can be helpful to work with the two-fold MB transforms by using integration by parts 
technique in the complex plane which is described in the previous section. Indeed, as we have 
derived in Eq. §§§ 



lim 

e 2 ->0,£i-s-0 



du dv x u y v M^ v) (e 1 ,e 2 ,e 3 ) 



c 



lim 

£ 2 ->-0,£i^0 



— / du dv x u y v D {u ^[l + e 3 ] + — y e2 I du dv x u y v D^[l - ej+ 

£l£2 Jc £ 2 £ 3 Jc 

—x £l I du dv x u y v D {u ' v) [l - e 2 ] = 
£i£3 



lim 

£2->0,£1^0 



c J 
du dv x u y v j\ y £2 [1 - £l ] + — J)M[l+£ 3 ] 

-£2^3 £l£2 
1 



+ x ei D {u - v) [l-e 2 ] 

£l£3 

From definition ([3]), we write the integrand explicitly 

— £>(»■") [1 + ea ] + —y E *D^[l - ei] + —x ei D {u ^[l - £l ] = 



£l£2 
T{-u)T(-v) 



£2£3 



^l£3 



1 r(ei + e 2 - u)r(ei + £ 2 - w)r 2 (l - e x - £2 + U + v) 
£l£2 



r(i + £ 1 + £ 2 )r(i-£ 1 -£ 2 ) 



1 



, T(£i - u)T(ei - v)T 2 (l - ei + u + v) 



£ 2 (£i + £ 2 ) a r(i + £l )r(i- £l ) 

1 r£l T(£ 2 - u)T{e 2 - v)T 2 (1 - £ 2 + u + v) 

r(i + £ 2 )r(i-£ 2 ) 



£l(£l + £2) 

The brackets in this formula can be re-written in the following form 



1 



£l£2 



rpi + £ 2 - M)r(£i + £ 2 - v)T 2 (l -e 1 -e 2 +u + v) 

r(i + £i + £ 2 )r(i-£i-£ 2 ) 

£2 r(ei - u)r( £l - v)T 2 {l - £l + u + v) 



r(i + £l )r(i- £l ) 



(17) 
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r(ei - u)r(ei - v)T 2 (l -ex+u + v) 



x ei T{-u)T(-v)T 2 (l + u + v) 



+ o(e 2 ) = 



r(l + ei)r(l-ei) 

—T(-u)r(-v)T 2 (l + u + v) [ip(-u) + i/>(-v) - 2ip{l + u + v) -lny] 

/ r(ei - u)T( £l - v)T 2 (l - £l +u + v^ (2) 

V r(i + ei)r(i- £l ) 

-lny (r( £l - u)r(ei - v)T 2 (l - £l +u + v))' ei=Q 

+—T(-u)T(-v)T 2 (l + u + v) [i/)(-u) + tp(-v) - 2^(1 + u + u) - lnx] 
£l 



El=0 



+ - 



1 (T{e 1 - u)Y(ei - v)Y 2 {\ -e 1 +u + v) 

r(i + ei)r(i-ei) 



(2) 



El=0 



-- In 2 xr(-u)T(-v)T 2 (l + u + v) + o(e 2 ) + o(ei). 

First of all, we demonstrate that all the poles in e% disappear. The corresponding contribution 
taking into account the factor r(— w)r(— v) has the form 

— / du dvx u y v T 2 {-u)T 2 {-v)T 2 {l + u + v) [{ip{-u) + ip(-v) - 2-0(1 + u + v)-\ny) 
£i Jc 

+ (tp(-u) + tp{-v) - 2^(1 + u + v) - lnx)} = 

— I du dvx u y v T 2 {-u)T 2 {-v)T 2 {\ + u + v) \(2ip(-u) - 2^(1 + u + v)-lnx) 
£i Jc 

+ {2ip(-v) - 2^(1 + u + v)- lny)] = 

-L I du dv(d u + d v ){x u y v T 2 (-u)T 2 (-v)T 2 (l + u + v)} = 0. 
£i Jc 



Thus, we have derived 



£2->0,l 



lim / du dv x u y v M^- v) (e 1 ,e 2 ,e 3 ) = 

>0,£l— >0 J (j 



du dv x u y v T{-u)T(-v) 



3 /T(£i - u)r(£i - v)T 2 (l -si+u + v) 
2 I ' 



(2) 



£1=0 



r(i + £i)r(i-£i) 

- lny (r(ei - u)T{e x - v)T 2 {\ -e x +u + w))^ =0 

~ in 2 .Tr(-M)r(--y)r 2 (i + u + v) 

As to the last two terms, we can consider the integral transformations 

du dv x 11 y v T(-u)T(-v) - lny (r(ei - u)r(ei - w)r 2 (l - e x + u + v))[ 

~ In 2 xT(-u)T(-v)T 2 (l + u + v) 

-lny j> dudv x u y v {ip(-u) + tp{-v) - 2^(1 + u + v)) r 2 (-u)r 2 (-v)r 2 (l + u + v) 



£1=0 
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-- In 2 x I dudv x u y v T 2 {~u)T 2 {-v)T 2 (l + u + v) = 
2 Jc 

~ lny / du dv x u y v (d u + d v ) T 2 {-u)T 2 {-v)T 2 {I + u + v) 
2 Jc 

~ In 2 x I dudv x u y v T 2 (-u)T 2 (-v)T 2 (l + u + v) = 
2 Jc 

— lnt/Oz + lny) * du dv x u y v T 2 (-u)T 2 (-v)T 2 (1 + u + v) 
2 Jc 

In 2 x I du dv x u y v T 2 (-u)r 2 (-v)T 2 {l + u + v). 
2 Jc 

Thus, the result for integral (IT71) is 

Km / du dv x u y v M^ v) (e u e 2 ,e 3 ) = 

5 £ " * y r( - u)r( ^ } [ r(i + gl) r(i- ei) ) £1=0 

-i (In 2 a; + In x In y + In 2 y) I du dv x u y v T 2 (-u)T 2 (-v)T 2 (l + u + v) = 
2 Jc 

| I dudv x u y v T(~u)T{-v) (T{e 1 - u)T( £l - v)T 2 (l -s 1 +u + v)) e =q 

-- (ln 2 x + lnxlniy + ln 2 t/) * du dv x u y v T 2 {-u)T 2 {-v)T 2 {l + u + v) 
2 Jc 

~ (T(l + ei)r(l - £i))^L I du dv x u y v T 2 (-u)T 2 (-v)T 2 (l + u + v). 

j c 

This formula can be developed further and be transformed to the form without derivatives of the 
Euler ip function in the integrand, 



du dv x u y v T(-u)T(-v) (Tfa - u)T{e x - v)T 2 {l - e x + u + v)y ' = 
c 61 



du dv x u y v [—d u ip(—u) — d v ip{—v) + d u ip{\ + U + v) 

'c 

+d v ip(l + u + v) + {ip{-u) + - 2^(1 + u + v)) 2 ] r 2 (-u)r 2 (-i;)r 2 (l + u + v) = 

du dv x u y v [(In x - 2ip(-u) + 2ip(l + u + v)) (ip{—u) - ip(l + u + v)) + 

< • 

+ (bay — 2ip(-v) + 2ip(l + u + v)) (tp(-v) -ip(l + u + v)) 

+(i>(-u) + tp(-v) - 2ij>(i + u + v)) 2 ] r 2 (-M)r 2 (-w)r 2 (i + u + v) = 

du dv x u y v [lnx {ip{—u) — ip(l + u + v)) + lny (ip(—v) — + u + v)) 

-(ip(-u) - ^(-v)) 2 ] r 2 (-u)r 2 (-u)r 2 (i + u + v) = 

hln 2 x + ln 2 y) <f du dv x u y v T 2 (-u)T 2 {-v)T 2 (l + u + v) 
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du dv x u y v {ip{-u) - i){-v)yT z (-u)T\-v)T z {l + u + v) 



The first term is proportional to the first UD function, and the result for integral (|17[) is 

lim f du dv x u y v M^ u ' v ' (£1, e 2 , £3) = 

£2— >0,£l— >0 J Q 

- In 2 - - ~ (r(l - e)r(l + e))fj ) / du dv x 11 y v T 2 (-u)T 2 {-v)T 2 (\ + u + v) - 



c 



du dv x u V v - ^{-v)f r 2 (-w)r 2 (-v)r 2 (i + u + v) 



(18) 



This result is related to formula (|16p via integration by parts in the complex planes of u and v. 



6.3 Explicit two-fold MB transform for the three-rung ladder 

We need to derive the same limit for M3 (ej.,£2,£3)- We show that the result is reduced to the 
Euler ij) function and its derivatives. According to Eq. ([±2)) , 



1™ . n M t' V \ £ l> £ 2' £ 3) 



£2— ►0,ei->0 



lim 

e 2 ->0,ei-M) 



-M, 



J 



(ei, £2, £3) t >' . 

£2£3 



-M. 



(u+si,v) 



lim 

e 2 — >0,ei->0 



£l£2 



-M. 



r(u+ei,v+e2) 



(£l,£2,£3) 



J 

£l£3 



-M. 



(u,v+e 2 ) 



£2(£l + £2) 



M 



r (u+£l,v) 



J 



£l(£l + £2) 



-M, 



(u,V+E2) 



fo) 



lim 



/ 1 



M. 



(u+ei ,'u+£ 2 ) 



(£i,£2,£3) - JAf^ 



£2->0,ei->0 \£i£2 

We use the notation 

(T(-e< - U )r(£, - v))f , (T(- Si - u)r(e> - v)) ( V 



M. 



(u+Ei,t)) 



(u,v) 



(19) 



(r(-e, - u)T( £l - v)) ( 



(4) 



for the corresponding derivatives with respect to £$ at the point £j = 0. We obtain for the second 
term of (|19p the expansion in terms of e± 



M ("+euv)^_ M (u,v) 



r 2 (i + u + v)' £ -± (r(- £l - U )r( 6l - «))< 3) r(-«)r(-«) 



£1 



(r(-£! - u)r( £l - v))' (r(-ei - u)r(ei - «)) 



(2) 



+ £l a^ (r(- £l - u)r( 6l - v))' n-u)r(-v) + -± (r(- ei - «)r( £l - «)) w r(-«)r(-«) 

-f (r(- £l - U )r( £l - v))f (r(- £l - u )r( 6l - «))(, 
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+4 2) £ ?(^ (r(- £l - U )r( ei - v))^ n-u)r(-v) 



\ (r(- £l - u )r( El - v))' (r(- £l - u)r( 6l - v))' ) + af > £ l±r 2 (- u )r*(-v) 



+ o(el). 



We use the notation 



a( £ ) = [r(i- £ )r(i + £ )]- 



» 



(a( £ ))^o 



Going back to the first term of Eq. (|f 91) we derive, by using Eq. ©, that 



1 r(-ei - u)r(-g 2 - f)r(e 2 - u)r( £l - iQr 2 (i + u + v) 

sie 2 r(l-e 3 )r(l + e 3 ) 

r(-ei - u)T(-v)T(-u)T(ei - v)T 2 {l + u + v) 



I 



S2(ei+e 2 ) r(f-ei)r(l + ei) 

1 r(-u)r(-e 2 - v)T{e 2 - u)T{-v)T 2 {\ + u + v) 



ei (ex + £ 2 ) 
A simple calculus gives that 



r(l-e 2 )r(l + e a ) 



M^ u+ei ' v+e2) ( £l ,e 2 ,e 3 ) - JM^ u+ei ' v) ( Sl ) 



e 2 r 2 (f + u + v) 



-_(r(- ei -u)r( El -«))^r(-«)r(-t;) 



+ o ( r (- £ i - u ) r ( £ i - u ))o ( r ( £ 2 ~ u)r(-e 2 - «)) 



(2) 



4 2) (r(- £l - U )r( ei - «))£ r(-«)r(-«) 



(r(- £l - «)r( El - v)i 4) r(-u)r(-v) 
+ £ j- (r(- £l - U )r( £l - «))( 2) (r( ea - u)r(- ea - w ))|, 2) 
+4 2) £l (i (r(- £l - u )r( E i - v))^ r(-u)r(-v) 



(r(- 



u)r( £l - «))' (r(- £l - «)r( ei - v))' ) + a< 4) £l ^r 2 (- u )r 2 (- v ) + ( £l ) 



o(e 2 ). 



Thus, we obtain the limit for M^ u ' v \si, s 2 , 63) when the e-terms are vanishing 



Jj m _ k n M 3 U '" ) ( e l' e 2' e 3) 
e 2 — >0,ei— >0 



lim 

£2^0, £1^0 \ £i£ 2 



1 



M. 



(u+£i,D+£ 2 ) 



(e 1 , S2 ,e 3 )-JM^' v \e 1 ) 



+ - 



1 r 



r 2 (i + u + u) 



i (r(- £l - u)r( £l - v))W (r(- £l - u)r( £l - 



(2) 
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~ (T(-ei - u)T(e 1 - „))< 3 > (r(-e x - u)r( 6l - 

+4 9) (-| ((r(-ei - t»)r( ei - «))o) 2 - r(-«)r(-«) (r(- £l - u)r( £l - v))^) 



(21) 



7 Conclusions 



In this paper we have explicitly calculated a new type of the two-fold MB integrals. The integrals 
of that type have been obtained due to recursive relations of the momentum loop integrals which 
are derived from the Belokurov-Usyukina loop reduction method. The recursion creates an infinite 
system of the integral relations, which allowed us to represent the MB transform of momentum 
(L + l)-loop integral as a linear combination of the MB transforms of momentum L-loop integrals. 
All the integrals have been reduced to the Appell's hypergeometric function corresponding to the 
one-loop momentum massless integral. In order to derive those relations, we needed to shift the 
complex variables of the MB transforms on the right hand sides of the integral relations. This shift 
is necessary in order to organize the same dependence on the parameters x and y in the integrands 
on the on the left and on the right hand sides. As the result, in the integral relations representation 
Q16p appears while representation (|18l) does not appear. 

At the first site, the results (fTgj) and (|2"Tj) are difficult to analyse, since on the right hand sides 
of them certain products of higher order derivatives of the Euler gamma functions 

-^n-e ~ u)T(e - v)\ E=0 (22) 

appear. These derivatives can be re-written as 
if 

— r(-e - u)r( e - v)\ £=0 = (flu - d v ) n T(-u)T(-v), (23) 

that makes them more convenient for use in the integration by parts. 

However, the integration by parts in the planes of complex variables u and v can serve to 
us just to prove the coincidence of two, at first view, different representations. The reason why 
the representation for MB transforms of UD functions should be a simple expansion in terms of 
derivatives (f2"2"| or (|23[) is in a recursive construction of the MB images. The recursive construction 
has nothing to do with quantum field theory or with theory of polylogarithms. We will prove this 
observation further. For thispurpose representation (|18l) is more useful than representation (| 16[) 
used in the integral relations [5 

As the first step, we consider the integral relation @ that can be found in Ref. [1] 

1 

2' 



J(l, 1, 1 - e) = ~(plY [x s J(l - e, 1 + e, 1) + f J(l + s.l-e, 1)] 



or, equivalently, 



j> du dv x u y v D iu ' v) [l,l,l-e] 



3 Actually, beyond the four rung in the ladder the difference between two types of representation disappears 
4 This relation can easily be derived via shifts in the planes of complex integration variables u and v in the MB 
transforms 
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- I dudv x u y v [x £ D^ u ' v) [1 - e, 1 + e, 1] + y £ D^ v) [1 + e, 1 - e, 1]) . 



(24) 



The benefits of representation (|24p are in more simple dependence on the variable e, namely, it 
presents in the numerator in two factors only, 



[l-e,l + E,l] = 



T{-u)T{-v)T{-£ - u)T(e - v)T 2 (l + u + v) 



r(l-e)r(l + e) 
For example, the equation for M 2 can be re-written as 



lim / du dv x u i/"M 2 M (ei,£2, £3) = 

>0,ei-yOJ c 



lim d> du dv x u y l 

£ 2 ^0,£i-i-0 Jq 



y s *D^[l - ei] + £>(»>") [1 + e 3 ] + 

£2£3 £l£2 £l£3 



B «i£>(«».«)[l_ e3 ] 



lim (/> dit x u y° 

£2— >0,£i — i-0 



a(£i) 

. £2£3 

a(e 3 ) 



y 62 [x £l T{-u - £i)r(-« + ex) + y ei T(-u + £i)T{-v - £l )} 



£i£2 
a ( £ a) rS1 
£l£3 " 



[x~ £3 r(-u + e 3 )r(-w - £3) + y- e3 T(~u - e 3 )r(-« + £3)] 
[x S2 T(-u - e 2 )T{-v + e 2 ) + y £2 T{-u + e 2 )T{-v - e 2 )] x 



xT(-u)T{-v)r 2 (l + u + v). 
In the limit of vanishing gj the value of this expression contains derivatives (j2"2"]l only, that is, 



£2-^0 



lim / du dv x u y v Mt' v) (£ 1 ,e 2 ,s 3 ) 

>0,£i^0/ c 



du dv x 11 y v T(~u)T{-v)T 2 {\ + u + v) 



^(a(e)r(-u-e)r(-v + e)) { Q 2) + 



1 In - (o(e)r(-u - e)T(-v + e))' + ~ In 2 -T(-u)T(-v) 

2 y 4 y 



This expression does not have a form of (|16l) or (|18p . however it can be transformed to that form 
by using the integration by parts in the complex planes of u and v. 
In integral relation (TT2"j) for M 3 , which is 

Mt v) (e 1 ,e 2 ,e 3 ) = 

1 -x-^y-^M^ v \ £l ,e 2 ,s 3 ) + J- x ~^ M^ v) (e,) + J- y -^M^\e 2 ), 



£l£2 



£2£3 



£l£3 



we substitute the structure found in (|24p. and again in the limit of vanishing we obtain an 
expansion in terms of derivatives (|22[) . that is, 

lim a> du dv x u y v Mi u ' v \si, e 2 , s 3 ) = 

£2->0,£l— >0 Jq 



du dv x u y v T(-u)T(-v)r 2 (l + u + v) 



^(a(e)r(-u-s)r(-v + s)i 4) + 
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I In - (o(e)r(-« - e)T(-v + e))f + \ In 2 - (a(e)r(-u - e)r(-« + e))f 
6 J/ 2 y 

' ' = ' " " ' (25) 



-^^(a(e)T(-u-e)T(-v + e))' 

By using the integration by part in the complex plane of u and v, we obtain coincidence with Eq. 

Following the recursive procedure for higher values of n, we will obtain the higher derivatives 
of type (|22[) multiplied by powers of ln^. At present, it is difficult to calculate the coefficients in 
front of the expansion terms. However, it is clear that the finiteness of the limit of vanishing £j 
has nothing to do with the poly logarithms, since instead of the Euler T function we can write any 
other smooth function in construction (|24[) and the limit remains finite. Thus, the infinite sum in 
the limit of vanishing Si of the quantities M^' v \ei, £2, £3) constructed from (|2"4")l has underlying 
integrable structure which can be uncovered by identifying coefficients in front of the terms in the 
expansion (|25[) for higher n. That integrable structure has nothing common with the MB transforms 
of polylogarithms and is based on the properties that can be studied by the basic methods of 
mathematical analysis. 
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